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The question of efficient multimode description of optical pulses is studied. We show that a rel- 
atively very small number of nonmonochromatic modes can be sufficient for a complete quantum 
description of pulses with Gaussian quadrature statistics. For example, a three-mode description 
was enough to reproduce the experimental data of photon number correlations in optical solitons 
[S. Spatter et al., Phys. Rev. Lett. 81, 786 (1998)]. This approach is very useful for a detailed 
understanding of squeezing properties of soliton pulses with the main potential for quantum com- 
munication with continuous variables. We show how homodyne detection and/or measurements of 
photon number correlations can be used to determine the quantum state of the multi-mode field. 
We also discuss a possible way of physical separation of the nonmonochromatic modes. 
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I. INTRODUCTION 

For a complete quantum description of an optical 
pulse, one has to use a multimode density matrix. The 
question is, how many modes are necessary for such a de- 
scription. If one works with monochromatic modes, then 
an infinite number of modes would be needed, which is 
impractical. On the other hand, one can construct differ- 
ent sets of modes as linear combinations of the monochro- 
matic modes 0. These modes (which are however non- 
monochromatic) may be more suitable for the quantum 
description of pulses. In particular, it would be useful to 
find such a modal structure so that the quantum state of 
the pulse can be described by the density operator g = 
Pcxc<X)Pvac, where g exc is some nontrivial density operator 
of a few modes and g V ac is the vacuum state density op- 
erator of all the remaining modes. Thus, we could work 
with a concise description of the pulse by means of (5 C xc- 

Our work is motivated by the question how to com- 
pletely describe the quantum state of solitons in optical 
fibers. Such soliton pulses are used in various schemes of 
quantum information processing (for a review see, e.g., 
||). Can a single- mode description of a soliton pulse 
be sufficient? When explaining squeezing in the Kerr 
medium, one sometimes plots a picture of phase space 
where a quantum uncertainty circle is deformed into an 
ellipse — a description which clearly corresponds to a sin- 
gle mode situation. This approach is useful to qualita- 
tively understand the generation of continuous variables 
entanglement as in ||, but do we use all the quantum 
features of a given pulse if we treat it as a single mode ob- 
ject? Should one, instead, work with many more modes 
of a pulse in the hope that each of them can become 
a useful resource for quantum communication? Multi- 
mode correlations of photon numbers in a soliton were 
observed H. To numerically calculate such correlations, 



quantum variables of a soliton were treated on a posi- 
tion grid of typically ~ 10 2 — 10 3 points ||, |(|. Are there 
hundreds of useful modes available in a pulse, or would 
a better choice of the mode functions show that only a 
few modes (perhaps just one?) are in some nontrivial 
quantum state? Or are just four quantum operators in- 
troduced by Haus and Lai Q sufficient to describe all the 
relevant phenomena? 

Knowing the answer to the question of what is the 
complete quantum description of a pulse would be very 
helpful in quantum information processing: one could 
fully utilize the squeezing and entanglement properties of 
our sources. After forming a soliton pulse in a fiber, one 
can determine its multimode quantum state in the most 
comprehensive way. One can then optimize the medium 
properties to achieve maximum squeezing, or maximum 
purity of an entangled state. Working with pairs of cor- 
related pulses, one could apply a proper measurement 
scheme and use entanglement criteria for multimode bi- 
partite Gaussian states || to check wheteher the pulse 
pair is entangled or separable. One can also better un- 
derstand the influence of the medium on the propagating 
pulse: provided that the output pulse is deformed, what 
happens with the quantum information carried by the 
pulse? Is it washed out by decoherence processes (or 
perhaps by an eavesdropper), or is it just unitarily trans- 
formed into other modes of the same pulse? A simple 
and correct measurement and description of the multi- 
mode state is highly desirable. 

This work is organized as follows. In Sec. [n] we intro- 
duce nonmonochromatic modes and deal with the trans- 
formations between different sets of quadrature opera- 
tors. Sec. Ill studies basic properties of the photon 



statistics of multimode fields: mean photon numbers and 
covariances between different modes. In Sec. IV we dis- 



cuss a homodyne scheme for a complete determination of 
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Gaussian states of nonmonochromatic multi-mode fields. 
In Sec. we compare the photon number squeezing 
available via applying a proper local oscillator modula- 
tion and via using spectral filtering. Sec. (v| suggests a 
way for an optimal selec tion of the nonmonochromatic 
mode functions. In Sec. VII we discuss a possible way 
how to physically separate individual nonmonochromatic 
modes of the pulse. Discussion and conclusion are pre- 
sented in Sec. [VIII . Most of the mathematical details are 
discussed in Appendixes [A|-[e], and in App. |] we compare 
our approach to that of Haus and Lai (?]] . 



II. NONMONOCHROMATIC MODES 

Under the term "mode" we understand a single de- 
gree of freedom of the electromagnetic field; a mode can 
be described by a pair of bosonic operators. It can 
be monochromatic (evolution described by a single fre- 
quency) or nonmonochromatic. Any state of the field 
can be treated in different mode decompositions. In a 
given decomposition the state is called single mode if all 
modes except of one have vacuum statistics of their op- 
erators, in the opposite case the state is multimode. In 
this section we deal with the transition between different 
mode decompositions. 



A. Bosonic operators 

Let us assume polarized optical field propagating in a 
given direction. The monochromatic modes of the field 
are denoted by the corresponding frequency to, and the 
bosonic operators of these modes satisfy the commuta- 
tion relations 

[a(u),tf(u/)] = 5((j-cj'), [fi(cj),fi(w')]=0. (1) 
Let us assume a complete orthonormal set of functions 



(2) 
(3) 



fm(^)f k (uj)duj = S mk , 

53 /jjH/fcCw) = 



We define a new set of operators b\ as 

h = J fk{oj)a(w)dw (4) 
which satisfy the commutation relations 

[h,bl}= 5 kk>, [b k ,b k ']=0. (5) 

Thus, the functions can be used to define a new set 

of nonmonochromatic modes. The inverse transformation 
of the nonmonochromatic modes to the monochromatic 
ones reads 

a(cu) = J2fk(")h. (6) 

k 



B. Example 

Let us consider a pulse with a normalized frequency 
envelope /(w); J | / (ui)\ 2 duj = 1. Let us define an or- 
thonormal system of mode functions as in Eqs. @ and 
(||) with fi(oj) = /(^)- Let the quantum state of the 
first mode &i be a coherent state = P\f3)i, 

and let all the other 6-modes be in vacuum 0)fe, bk\0)k — 
0, k > 1. By means of the transformation (p) we find that 
using the monochromatic modes, the quantum state is a 
multi-mode coherent state rj w Even though 

the single-mode and multimode coherent states are re- 
lated to the same object, i.e., a pulse in a coherent state, 
the single-mode description is clearly more convenient 
for handling and for understanding the field structure. 
Of course, for other than coherent states the transfor- 
mations of the state in different mode-systems are not 
so straightforward, but still can provide us with a very 
convenient way of the quantum state description. 



C. Quadrature operators 

We define the hermitian quadrature operators x(ui), 
p(oj) as 



x(uj) = (a(u) +a f (o;)) , 
PM = —7= (a{u)-a\u)) 



(7) 

(8) 



and similarly the 6-mode quadrature operators X k , P k as 
Xk = -j=(b k + bl), (9) 

^ = ^(V*£). do) 

These operators obey the same commutation relations as 
the quantum mechanical position and momentum oper- 
ators with h = 1, i.e., 



[x(u),p(u')) =iS(u-u}'), 
[x(u),x(uj')) = ]p(u),p{u)')] = 0, 



and 



X k ,P k 



i5 kk >, 



Afc, X k > 



Pk,P k - 



0. 



(11) 
(12) 



(13) 
(14) 



The transformation between the two sets of quadrature 
operators can be written in the form 



H=X,P k 



(15) 
(16) 
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with £ = x,p, and /z 
matrix Z is 



X, P, and the transformation 



Re / fc (u,) (17) 

z*>) = -z£>) = im (18) 

As can be checked, matrix Z satisfies the orthogonality 
relations 



E 



(19) 
(20) 



Sometimes it is useful to work with a discrete set of fre- 
quencies u>j corresponding to frequency bins of width Alu. 
Quadrature £j of the jth bin is obtained as 



fj,=X,P k 



(21) 



(22) 



For brevity, we can join the quadratures x(w) and p(u>) 
into a single vector £, and similarly the quadratures X k 
and Pfc into a single vector (a, so that the transforma- 
tions jlq ) and ( [Tq ) [or (|2l|)l are written in the matrix 
multiplication form 



i = z T p. 



(23) 



Here the superscript T means matrix transposition, and 
the rows of the matrix Z are represented by indices /j, 
and k, and the columns by the indices £ and to (or j). 



D. Mean quadratures and variances 

Important properties of multimode quantum states are 
given by the vectors of the quadrature mean values, £j = 

Tt[££j] and ~p k = Tr[gp, k ], and by the variance matrices 
V and V, given as 

Wfc' = 5({(w-Mit),K'-A0})» ( 24 ) 

= 5<{fe-?j).(^-^)}> (25) 

where {. . .} stands for the anticommutator, {A, i?} = 
AS + BA, and (. . .) represents averaging (A) = Tr(gA). 
The relationship between these quantities can be written 
in the matrix multiplication form as 



I = z T -p, 

V = ZV'Z T , 
V' = Z T VZ. 



(26) 
(27) 
(28) 
(29) 



In general, to transform between ~p, V, and £, V', one 
has to know the complete transformation matrix Z , i.e., 
the whole set of mode functions However, if we 

assume that only a few nonmonochromatic modes are 
excited while the rest is in vacuum state, then the explicit 
form of the empty mode functions does n ot pl ay any role. 
(It was also illustrated in the Example II B where only 
one mode function f(to) was enough to transform the 
state into the multi-mode case.) This can be used for 
a simplified calculation of the transformed quantities, as 
discussed in App. [A|. 



III. PHOTON STATISTICS OF MULTIMODE 
GAUSSIAN STATES 

Among all possible states, the class of Gaussian states 
is one of the most important — both from the theoreti- 
cal and experimental point of view. The basic property 
of these states is that their Wigner function (as well as 
the Q-function or the characteristic function) have Gaus- 
sian form (see, e.g., & |[ |l(J). Examples of Gaussian 
states are coherent states, squeezed coherent states, ther- 
mal states, and squeezed thermal states. Gaussian states 
are typically observed in most experiments with optical 
pulses. In this paper we confine ourselves to Gaussian 
states; the main advantage used here is that a Gaussian 
state remains Gaussian in any mode decomposition and 
a relatively small number of parameters is necessary for 
its description. 



A. Mean photon numbers and variances 

Let us first study the relation between quadrature mo- 
ments and photon number moments in discrete modes. 
The photon number operator in kth mode can be ex- 
pressed by means of the quadrature operators as 



1 



■pI-i) 



(30) 



(31) 



The mean photon number in kth mode is thus 

(n k ) = \ «4> + (Pi) 1) , 
and the mean product of photon numbers is 

(n k m) = i {{ilii) + {x 2 k pi) + {plpj) + (p 2 k xf) 

- (4) - (Pi) - - (Pt) + 1) • (32) 

As discussed in App. |B| for Gaussian states all the 
quadrature moments on right hand sides of Eqs. ( |3l| ) 
( p2j ) can be expressed usin g the quadrature means and 
variances as m Eqs. ® - flB10|) , The resulting photon 
number moments can be used to calculate the photon 
number covariances 



cov (n k ,ni) = (n k ni) - (n k )(m), 



(33) 
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which can be written in terms of the quadrature means 
and variances as 

cov (n k , m) = x k xiV.^. k xl + x k yiV' xk pl 



4-Z (V a A- V' 2 -4- V' 2 + V' 2 \ 
9 V xk,xl ' v pk.xl ' v xk.pl > v pk,plJ 



04) 



with the special case An^ = cov(n k ,n k ). A very im- 
portant role is played by the so called normally ordered 
covariancc 

C k i = (■ An fc An ; :) = cov(n k ,ni) - 5 kl (n k ). (35) 

The normally ordered covariance is used to define the 
normalized correlation matrix as in 101 , 



C, 



(n) 
kl 



An k Ani 



(36) 



Values of the correlation matrix measured in Q are 
shown in Fig. ^a. In the limit of continuous frequency 
modes, one can define photon number density (n(w)), 
variance density (An 2 (w)) = (n(u>)), and the normally 
ordered covariance C(u>,u>') with its normalized version 
(/'"'(UjW'). Since a Gaussian state remains Gaussian in 
any mode decomposition, and since the parameters of a 
Gaussian state £ and V can easily be transformed from 
one mode decomposition into another, one can also cal- 
culate the photon number mean values and correlations 
in arbitrary mode decomposition. 



B. Photon correlations in multimode pulses 




Wavelength ( nm ) 
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FIG. 1: Photon correlation matrix as a function of 

the wavelengths of the pulse spectrum. Results measured in 
[jl| with omitted elements with very large fluctuations (> 0.8) 
(a), and the best fit using a reco nstructed three-mode variance 
matrix Vxk,xk' ( see Sec. HID) (6). 



If we try to reproduce the experimental results of j|] 
as truly as possible with as few modes as possible, the 
first attempt would be to start with just a single non- 
monochromatic mode. However, as shown in Appendix 
|(| the covariance matrix cov(n.fc, ni) of such a field would 
have the same sign for all values k and I. The sign would 
be positive if the single-mode state is super-Poissonian, or 
negative if the single-mode state is sub-Poissonian. On 
the other hand, the covariance matrix measured in Q 
contains both positive and negative elements. Therefore, 
the quantum state of the soliton pulses measured in Q 
cannot be described as a single nonmonochromatic mode. 

Exact analytical expressions are rather lengthy if more 
than one mode are excited. However, relatively simple re- 
lations can be obtained if the following requirements are 
met: (i) The coherent amplitude of at least one of the 
modes is much greater than any of the variance matrix el- 
ements, (ii) only the X-quadratures have non-zero mean 
values, i.e., P k — for all k, and (Hi) the mode functions 
are real. While the condition (i) is generally valid for all 
strong pulses, the conditions (ii) and (Hi) are related to 
our choice of mode functions and their generalization is 
straightforward. In Eq. (fM) the terms containing p k ,i 



vanish and the first term on the right is dominant; thus 
Eq. (|34|) becomes 



cov (n k ,ni) w x k xiVL xl . 



(37) 



To calculate this quantity using the nonmonochromatic 
modes, we write 



x k 



= X] .fn(u k )X n \/~KL 



n<N 



(38) 



(see Eqs. (|27[), fllTj), and (|22|)), and for V^ k xl we use Eq. 
( A7 ) from Appendix The mean photon number is 



(n k ) ~ - ^ f m (ui k )f n (uj k )X m X n Auj, 



(39) 



using Eqs. (|3l|), (jBq), and the conditions (i)~(iii). The 
normally ordered covariance ( |35| ) thus becomes 

Sr, 



m,n<iV 
X 



C k l — X fm(u k )fn(ui) \ Vxm,Xn 7p" 

fm>(Uk)fn'(ui)X m .X n ,Acj 2 . (40) 



m'n'<N 
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Since in very narrow frequency bins Alu, the mean photon 
number and photon number variance are approximately 
equal, (An|) ~ (rife) <C 1, Eq. ( |39| ) can be used also 
to obtain (An?), and the normalized correlation matrix 
(pq) becomes 



fm(Uk)fn(u>l) (v Xm ,Xn - ^ ) A.M J I ! 
m,n<N ^ 



Taking Acj — ► 0, we can work with the continuous quan- 
tity 



C (n) (^V) : 



= 2 £ 

m,n<N 



'Xm,Xn 



~2~ 



(42) 



In Eq. ( |42| ) only the X-elements of the quadrature vari- 
ances V occur. In other words, the measured photon 
number covariances are only influenced by the covari- 
ances of quadratures which are in phase with the quadra- 
ture of the strongly excited modes. 



C. Reconstruction of the quadrature variances 
from the photon number covariances 

In the preceding subsection we have seen how to calcu- 
late the normalized photon covariance from the quadra- 
ture variance matrix of the nonmonochromatic modes. 
We can also consider an inverse problem. Let us assume 
that the mode functions /fe(k>) are known. Let us also 
assume that the normally ordered covariances C(lu,lu') 
are measured as in Q so that the normalized covariance 
(lu, u') can be determined. As can be seen from Eq. 
( [42| ) , this quantity is a linear combination of the products 
of the mode functions /fc(w). The coefficients in the lin- 
ear combination are elements of the quadrature variance 
V. Thus, to obtain the elements of V, one can use the 
orthonormality of the mode functions /fe(u;) to invert the 
linear dependence. We find 



Vxk. 



1 



Xk> 



>kk' 



(lu, u/)f k (w)fc (w'JdwtL/. (43) 



Thus, from the spectral covariances of the photon num- 
bers one can reconstruct the N(N + l)/2 elements of the 
quadrature variance matrix, out of the total number of 
N(2N + 1) independent elements. 

Since Eq. ( [43] ) is linear, it allows for a direct esti- 
mation of the reconstruction error (see, e.g., 0). If 
the normalized covariances were measured with preci- 
sion AC^(lu, lu'), then the error in reconstructing the 
elements of V can be estimated as 



AC {n)2 (uu,uu') 



xfl(uu)f 2 k ,(uu')dLudLu'. 



(44) 



Note that here, for the sake of brevity, we have assumed 
uncorrelated errors in the elements of (lu, lu'); deriva- 
tion of a more general formula taking into account cor- 
related errors is straightforward. 

Let us stress that the formulas derived in here and in 
Sec. Ill B are valid for very narrow frequency bins where 
(An 2 ,) 5 (rife) -C 1. However, in real experiments, it 
is often necessary to work with wider bins in which the 
photon statistics can differ from Poissonian so that using 
Eqs. ( |42] ) and (fl3|) could cause significant errors. This 
situation can be easily taken into account by substituting 
the proper expression for (An 2 .) in the formulas connect- 
ing the photon correlation matrix with the quadrature 
variance matrix. The equations (42) and J43| ) then be- 
come slightly more involved; since this generalization is 
straightforward, we do not include the formulas in this 
text. 



D. Application to experimental data 

To illustrate our method, we have used the experimen- 
tal data obtained in (reproduced in Fig. |l]a). They 
were measured for soliton pulses propagating in a 2.7 m 
optical fiber. We have used a preliminary set of mode 
functions (starting with a sech function as the basic shape 
of the soliton) and applied the reconstruction formula 
(p3|) . We have diagonalized the resulting variance ma- 
trix Vxk,xk' and found that only three eigenvalues are 
substantially different from the vacuum value 1/2. The 
corresponding eigenvectors can be used to construct a 
new set of mode functions (see Fig. |^). In this set, the 
X-quadratures are independent of each other. Thus, we 
have found that with respect to the quadratures which 
are in phase with the coherent amplitude (X-quadratures 
in our convention), the pulse measured in Q can effec- 
tively be described as a three mode field, the quadra- 
ture variances being Vxi,xi ~ 0.29, Vx2,X2 ~ 1-39, 
and Vx3,x3 ~ 2.69 with the off-diagonal elements being 
zero. The smallest eigenvalue corresponds to squeezing 
—2.35 dB. The reconstructed variance matrix Vxk.xk' 
can be used to calculate back the photon number corre- 
lation matrix Cjy ; see Fig. 0b. 

Since the photon number correlations were measured 
with limited precision, these results suffer from errors. 
To estimate the precision of our results, we have Monte- 
Carlo generated 1000 "experimental" matrices c|; with 
elements fluctuating with errors given by the error es- 
timates of the original experiment. The squeezing value 
then fluctuated between —2.2 dB and —4.1 dB, with most 
results centered around —3.1 dB. Even though the error 
is too large to make a definite statement, these result 
suggest that larger squeezing is available than the mea- 
sured —2.5 dB of this setup with spectral filtering. 
Let us note that the shape of the mode function corre- 
sponding to the squeezed quadrature (full line in Fig. |^) 
resembles the spectral filtering approach: the contribu- 
tion from the middle of the spectrum is enhanced while 
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FIG. 2: Mode functions for which the Vxk,Xk' matrix of the 
pulse measured in Q is diagonal. The insets show the fluctu- 
ations of the X-quadratures in comparison with the vacuum 
fluctuation (dashed line). In the main figure, the full line cor- 
responds to the mode function f\ with the squeezed quadra- 
ture, Vxi,xi ~ 0.29, dashed line to /2 with Vx2,x2 ~ 1.39, 
and the dash-dotted line to fs with Vxz,x3 ~ 2.69. 



the outer parts are suppressed. The question of photon 
number squeezing availability via spectral filtering and 
via local oscillator functions is studied in more detail in 
Sec. 0. 



E. Conclusion 

We have seen that one mode description of a pulse 
is not sufficient to explain the experimental data of dJ, 
whereas a three mode description gives a good agree- 
ment with the measurements. However, we cannot con- 
clude from this that the pulse does not contain more 
than three modes. In other modes the P-quadratures 
can be excited which do not influence the observed pho- 
ton statistics. Also, the measurement noise was too high 
so that weak excitation of some additional modes might 
be undistinguished from the data noise backgraund. 

It is interesting to compare our approach with that of 
Haus and Lai (sec also |ll], |lj, ) . In their case the 
quantum field is decomposed into a "soliton" part and a 
"continuum" part. The soliton field is described by four 
operators Art, A9, Ax, and Ap, related to the soliton 
energy, phase, position, and velocity, respectively. As 
shown in App. IF], these operators can be expressed by 
our quadrature operators of four modes, provided that 
the mode functions are properly selected. However, only 
two of these operators (An and Ax) related to two our 
modes (/i and $2) influence the photon statistics. Thus, 
confining ourselves to the four soliton operators of |7j 
would not be sufficient to describe the observed phenom- 
ena. To apply the formalism of [jfl, one would have to 



work with the full set of the soliton and continuum op- 
erators. 



IV. COMPLETE DETERMINATION OF THE 
VARIANCE MATRIX 

Even though one can obtain full information about the 
X-quadratures from the spectral correlations of photon 
numbers, one has no access to the variances Vpk.pw and 
Vxk.Pk' ■ To get full information on the multimode Gaus- 
sian quantum state, one has to perform phase dependent 
measurements. Homodyne detection is one example of 
such a measurement; recently it was studied how to ap- 
ply homodyne detection scheme to the quantum state 
reconstruction of a multimode optical field (lj, [l5|. Be- 
cause we confine ourselves to the Gaussian states, the 
task is easier than reconstruction of a general quantum 
state. 

To find the variance matrix V, one can use the scheme 
of [|L5j| with the local oscillator pulses shaped to the form 
of weighted combinations of the mode functions. It is 
very useful to subtract the coherent amplitude of the 
pulse by using a balanced Sagnac interferometer |lq ] as in 
Fig. [|. Here, two counterpropagating identical squeezed 
pulses are interfering at a 50%/50% beam splitter. In 
one of the outputs of the beam splitter the pulses in- 
terfere constructively and form a bright pulse, which 
is then used to form the local oscillator. In the other 
output where the pulses interfere destructively, squeezed 
vacuum (or a more general field with zero mean ampli- 
tude) is formed. The squeezed vacuum pulse has the 
same quadrature variance matrix V as each of the two 
counterpropagating pulses. 

The bright pulse is shaped interferometrically so that 
its envelope has the form of different combinations of the 
mode functions Such a pulse is then used as a 

local oscillator in a balanced homodyne detector. Thus, 
if the local oscillator is /fe(w), one can obtain the value 
of Vxk,xk, whereas with the local oscillator ifk(u) one 
can obtain Vpk,Pk- Knowing these values and using the 
local oscillator pulse of the form /fc(w) + /fe'(w) one can 
obtain the value of Vxk.xk', by using the local oscilla- 
tor form i(fk(u>) + fk>(u>)) one can obtain the value of 
Vpk,pk', and by using the form /fc(w) + ifk>(uj) one can 
obtain the value of Vxm.Pm 1 ■ Altogether, ^(2^ + 1) dif- 
ferent forms of the local oscillator are sufficient to obtain 
all the N(2N + 1) independent elements of the variance 
matrix V. Let us note that to increase the precision of 
the measurements, one can increase the number of differ- 
ent phases of the local oscillator; a Maximum-Likelihood 
method for parameter estimations of a single mode Gaus- 
sian state using many phases has been discussed in Jl7[ . 

Having found the variance matrix, one has all the infor- 
mation about the Gaussian state. Then we can immedi- 
ately see, e.g., what is the maximum available squeezing 
of the state: it is the value corresponding to the minimum 
eigenvalue of V. Typically, this value will be smaller than 
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FIG. 3: Scheme to measure all the elements of the quadra- 
ture variance matrix. The Sagnac interferometer produces in 
one output squeezed vacuum and in the other output a bright 
pulse corresponding to one of the mode functions (say fi(ui)). 
The pulse shaping device transforms the bright pulse enve- 
lope into different combinations of the mode functions. Using 
the resulting bright pulse as a local oscillator in a balanced 
homodyne detector, one can reconstruct the quadrature vari- 
ance matrix V . 



the minimum diagonal element of V, which means that 
the optimum squeezing is shared among different modes. 
For example, there is a quadrature squeezing in the fun- 
damental mode /i due to the Kerr effect (assuming mode 
functions as in App. |]), but the X\ quadrature is also 
correlated to the X 3 quadrature, because of the corre- 
lation between the pulse width and energy. Therefore, 
better squeezing can be expected to occur in a combi- 
nation of the modes 1 and 3 than in the isolated mode 
1. Thus, if the scheme is used to produce squeezing, the 
measurement scheme can serve as a tool for a selection 
of the optimum local oscillator. 

One can also see how much are individual modes en- 
tangled with each other. To take full advantage of this 
knowledge one has to be able to separate individual 
modes from each other and distribute them among Alice, 
Bob, and other entanglement consu mers. This is, how- 
ever, a rather nontrivial task; in Sec. VII we will briefly 
mention a possible approach to its solution. 



V. PHOTON NUMBER SQUEEZING VIA 
LOCAL OSCILLATOR MODULATION VS. 
SPECTRAL FILTERING 



Let us assume that we want to prepare a pulse with 
the maximum photon number squeezing, i.e., the pho- 
ton number fluctuates as little as possible. To quantify 
the photon number squeezing, one uses the Mandel Q- 



parameter |18 defined as 

_ (An 2 ) - (n) 



(45) 



where n refers to the photon number of the entire pulse. 
This quantity is negative for sub-Poissonian field and pos- 
itive for super-Poissonian fields. It was suggested in sev- 
eral works Jl9| that spectral filtering of the pulse can 
lead to improvement of the photon number squeezing by 
blocking frequency bands with correlated photon fluctu- 
ations and letting through frequency bands with anti- 
correlated photon numbers. Here we show that a proper 
modulation of the coherent amplitude (playing the role of 
the local oscillator) leads to the optimum squeezing which 
cannot be overcome by the spectral filtering method. 



A. Local oscillator modulation 

Let us first assume that the quadrature variance ma- 
trix V is fixed while we can modulate the mean values 
of the quadratures X n and P n ; this corresponds to the 
experimental situation as in Sec. IV and Fig. ^. The 
photon number variance is 



(An 2 ) =^cov {n k ,ni). 



(46) 



k,l 



Assuming as in Sec. [II B that the coherent amplitudes 
are much greater than the variance matrix elements, and 
that the mode functions are real (generalization to com- 
plex functions being straightforward), we can express the 
photon number variance as 

(An 2 ) = £ (XmX n Vx ■m.Xn 

m,n<N 

+ 2XmPnVxm,Pn + PjnPnVPm.Pn) ■ (47) 



The mean photon number is 

{»> = £<»*> = 5 £ {xl+Pi 



2 ) 

n I > 



(48) 



so that the Q-parameter of Eq. ( pf5[ ) can be written, after 
some algebra, as 



Q = 2 £ 

m,n<N 



,Xr 



Xn 



m.P 



Pn 



where 



X,, 



X,„ 



J2n<N ip^n + Pn) 
p 

^2n<N (-^n + Pn) 



(49) 

(50) 
(51) 



8 



Thus the Q~P arame ter can be expressed in the matrix 
multiplication form as 



Q = 2^X T ,P T ) 



Vxx Vxp \ ( X 



V PX Vpp \ P 



- L (52) 



where X and P are column vectors with the X m and P m 
elements, and Vxx, Vxp, and Vpp are the correspond- 
ing submatrices of the variance matrix V. The variance 
matrix is multiplied from the left and from the right with 
a unit vector, so that the Q-parameter is limited by 



2 | ^ (min) - - J < Q < 2 (v (r 



(53) 



where y( mm ) is the minimum eigenvalue of the variance 
matrix V and l/( max ) is its maximum eigenvalue. The 
minimum value of Q is reached if the vector of mean 
values of X and P is the eigenvector of V corresponding 
to its minimum eigenvalue. Setting the components of 
this vector is possible using the scheme of Fig. ||. 



B. Spectral filtering 

Let us assume a spectral filtering function < c(lu) < 
1: the frequency component is completely transmitted if 
c{lo) — 1 and it is completely blocked if c(u>) = 0. This 
function transforms the mean quadratures and variances 
as 



PI 



^2n<N c mnPn 



with 



C 2 (u>)f n (u))f n > {U>)du>. 



(61) 



(62) 



Again, the Q-parameter is calculated as a matrix product 
of the variance V multiplied from the left and from the 
right by a vector. This time, however, the vector is not 
of the unit length, so that Q* is limited by 

2 A 2 (v {min) - < Q f < 2 A 2 jV max > - ij , (63) 

where A 2 is the square of the magnitude of the multiply- 
ing vector, 



a 2 =y.(H 2 + p1 2 



i<N 



(64) 



To show that Q' can never be smaller than the mini- 
mum value achievable by local oscillator modulation, it 
is enough to show that A 2 < 1, i.e., that the magnitude 
of the multiplying vector is not bigger than one. Proof 
of this inequality is shown in Appendix 0. 



Conclusion 



-f 

Pk 



v'( 

pk,pl 



c{oj k )p k , 



Ski 



(54) 
(55) 

(56) 



ck)^ M c( Wi ) + y[l~C 2 K)] 

cK)V^, pi cM + ^ [1 - c 2 K)] , (57) 



V 



If 



xk,pl 



cMV^ p1 c{uji) 



(58) 



The Sm terms in Eqs. (|56j) and ( |57| ) follow from the 
quantum mechanical nature of the quadratures: partially 
blocking a frequency component means that the corre- 
sponding field is mixed with vacuum. 

Eqs. (|54|)-([58l) can be used to determine (An 2 ) and 
(n) as in the preceding subsection, and thus to find the 
Q-parameter of the filtered field. After a straightforward 
algebra, one can express the new Q-parameter as 



Q f = * E 



i<N 



xLivx 



m.Xn 



+2XlV Xm .PrM + PI Vp 



2 



where 



X s 



"yir, n'<N Cnn f {X n X n * ~\~ P n P n > 



(59) 



(60) 



We can see that no spectral filtering can improve the 
photon number squeezing below the minimum eigenvalue 
of the quadrature variance matrix, which is available via 
the local oscillator modulation approach. Let us stress 
that our method of finding the optimum local oscillator 
function is very simple and straightforward since it is 
based on linear algebra. It was suggested recently to use 
an adaptive algorithm to optimize the pulse shape for 
achieving optimum photon number squeezing pp| . In this 
method, the optimum was reached after 20 000 iterations. 
In our case, provided that the pulse is sufficiently well 
described by N = 5 modes, A^(2iV+ 1) = 55 iterations is 
enough. 



VI. SELECTION OF APPROPRIATE MODE 
FUNCTIONS 

So far it was assumed that the set of mode functions 
was given and the question was about the statistics of the 
corresponding quadratures. But how should these mode 
functions be selected? 

In principle, any orthogonal set of mode functions can 
be used for description of the pulse statistics. However, 
since the aim is to reduce the number of quantum vari- 
ables necessary to describe the pulse, the functions should 
be carefully chosen. One possibility is to start from the- 
ory and assume some particular shape of the pulse - e.g., 
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a hyperbolic secant soliton and to construct the orthog- 
onal set from the typical perturbations. An example is 
given in App. ||], Eqs. (F15)-( [fTs| ) where relationship to 
the soliton quantum fluctuations approach by Haus and 
Lai Q is discussed. 

Another approach does not assume any particular 
mode function form and is related directly to the ex- 
periment. It is useful to assume that only one mode has 
a nonzero coherent amplitude and a very small number 
of modes have quadrature fluctuations substantially dif- 
ferent from the vacuum values. To find the optimum set 
of mode functions, one can use the following procedure. 

1. Select fi(ui) as the measured classical envelope of 
the pulse (see App. [I]). 

2. Determine the minimum value of variance to be still 
considered as different from vacuum. 

3. Construct a temporary set of orthogonal functions 
/f mp) , k = 1 . . . N' with f[ temp) = fx chosen in |. 
The size N' of the temporary set should reasonably 
correspond to the experimental conditions. 

4. Perform the measurement of the 2N' x 2N' vari- 
ance matrix with the temporary set of mode 
functions. 

5. Construct the reduced matrix V 1 - 2 ^ from by ex- 
cluding rows and columns referring to quadratures 
of mode 1. Let indexes 1 . . .N' — 1 correspond to 
the X quadratures, and indexes N' ... N' — 2 cor- 
respond to the P quadratures. 

6. Diagonalize V (2 ^ to get = WV (2 ^W T with 
W an orthogonal matrix. Let us choose W such 
that V^p is the largest element of V^ 2 



£ w w lk WuV h 



(2) 
k,l ■ 



v n — 



7. The function f 2 is constructed as 



h = 



N'-l 

E 



{W lk +iW lk+N ,. x )f, 



(temp) 

(fc+i) • 



(65) 



It can be checked that the variance of the X 
quadrature corresponding to this mode function is 



lH2) 
v ll ■ 



A new temporary set of N' — 2 mode functions is 
constructed from the old one as an orthogonal com- 
plement to f\ and f%. By means of the transfor- 
mation connecting the new temporary set of mode 
functions to the initial one, calculate the corre- 
sponding variance matrix. Construct reduced ma- 
trix by excluding rows and columns referring 
to quadratures of modes 1 and 2. Diagonalize V^ 3 -* 
and find mode function f$ in the same way as in |6| 
and 0. Note that the variances of the quadratures 
of mode 3 are smaller than the variance (X 2 ). 



9. Repeat this procedure of redefining the temporary 
mode set, transforming and diagonalizing the vari- 
ance matrix, and defining a new mode function. 
After each repetition, the maximum variance of the 
(n + l)st mode is smaller than the maximum vari- 
ance of the nth mode. If for some n = N the vari- 
ance of mode n + 1 is sufficiently close to 1/2 (as de- 
fined in item ||), the quantum state of the (n + l)st 
mode is indistinguishable from vacuum. The pulse 
can be described with the given precision as an N- 
mode object. 

Let us note that the procedure of redefining the tem- 
porary mode set and rediagonalizing the reduced vari- 
ance matrix each time when a new mode function is con- 
structed is necessary. It is not possible, as one might be 
tempted, to find a set of mode functions simply by diago- 
nalizing the measured variance matrix, since the SO(2A^) 
transformation corresponding to diagonalization is gener- 
ally not a canonical transformation (see, e.g. || for more 
details). On the other hand, in the special case of pure 
Gaussian states which are specified by N(N + 1) real pa- 
rameters one can find a set of uncorrelated modes. In this 
case our procedure would be terminated after the first 
diagonalization. The diagonalization of pure Gaussian 
multimode states into uncorrelated modes has recently 
been studied in plf . 

The procedure of operational construction of mode 
functions as described above is, of course, not the only 
possible. It is, however, very useful for finding the mini- 
mum subspace of mode functions sufficient for the pulse 
description. Other sets of mode functions can be selected 
such that the variance matrix takes some special shape, 
e.g., some of the "canonical" forms studied in B. 



VII. SEPARATION OF MODES 

It may be very useful to separate individual non- 
monochromatic modes which form the pulse. The pos- 
sibility of obtaining nonclassical correlations of optical 
pulses by partitioning the pulse in the spectral region 
was suggested in |p2| . Similarly as with spectral filter- 
ing, this approach is not necessarily the optimum one for 
obtaining maximum entanglement from the source. 

The basic idea for obtaining the optimum separation 
is to send different (nonmonochromatic) modes into dif- 
ferent channels by using a special unitary transforma- 
tion among them. It was shown in p3~f that any dis- 
crete unitary operator can be constructed interferomet- 
rically. For this purpose we suggest to apply a scheme 
as in Fig. [|. In the first step one decomposes the pulse 
into quasimonochromatic components. The frequencies 
of these components are then shifted by means of acusto- 
optical modulators (AOMs) so that each channel has the 
same central frequency. The channels then interfere on 
a 2iV-port consisting of beam splitters and mirrors. By 
a proper choice of the 27V-port parameters one can man- 
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age that most of a pulse in the /fc(w) mode leaves the 
27V-port in the kth output channel. 




FIG. 4: Scheme to separate individual nonmonochromatic 
modes. The pulse is first decomposed on a grating into quasi- 
monochromatic channels. The frequency of each channel is 
then shifted by means of an AOM to the central frequency 
of the pulse. The resulting modes (with the same central 
frequency) then interfere on a 2iV-port. 

One can also optimize the 2iV-port parameters to pre- 
pare a set of optimally entangled K modes (multipartite 
entanglement), sent to different channels. This would be 
a generalization of the proposal of |22| for partitioning 
soliton pulses to generate entangled states. 

Let us note that the spectral filtering of solitons used 
to produce photon number squeezing p9| is a special kind 
of mode separation. However, as shown in Sec. in this 
case the separation is not optimized with respect to all 
relevant degrees of freedom. Therefore, better results can 
be expected in our general approach. 



VIII. DISCUSSION AND CONCLUSION 

To summarize our results we can state that: 

(i) For a complete description of the experimental re- 
sults measured on optical solitons it appears that a very 
small number of nonmonochromatic modes is enough. 
The results cannot be described by means of a single- 
mode field, but already three nonmonochromatic modes 
were sufficient to reproduce the experimental data of [Q . 
With our approach it is easy to interpret the "butterfly" 
pattern of the measured photon covariances: three dif- 
ferent non-monochromatic modes overlap and contribute 
with their fluctuating in-phase quadratures to the photon 
statistics (see Fig. ||). 

(m) Provided that the pulses are Gaussian (which 
seems to be a relevant assumption for most of the exper- 
imental situations), a complete quantum description of 
the pulse can be done by means of a multimode variance 
matrix. The elements of this matrix can be determined 
by means of homodyne detection with specially shaped 
local oscillator pulses. 

(Hi) Concrete form of the mode functions is a matter 
of choice. If the aim is to find the smallest number of 



modes, an operational method for constructing the mode 
functions is provided. In this case the first mode contains 
the coherent amplitude, whereas the rest of the modes 
have zero mean fields and their field variances decrease 
with increasing mode index. By selecting the minimum 
set of modes, one can substantially reduce the number 
of parameters necessary for a complete description of the 
physical situation. 

(iv) Knowledge of the mode structure of the pulse and 
of the corresponding quantum state will be very useful 
for quantum information purposes: one can select the op- 
timum shape of the local oscillator pulse to detect max- 
imum squeezing or entanglement. This optimum finding 
is very straightforward and much faster than adaptive al- 
gorithms as in J2(|. One can also relatively easily study 
the influence of the medium on the propagated pulses 
and on the quantum information they carry. By mea- 
surement of the multimode quantum state of the input 
and output pulses we can find the von Neumann entropy 
of the states. This would enable us to tell whether the 
observed pulse deformation corresponds most probably 
to a unitary evolution or rather to decay and dephasing, 
possibly caused by an eavesdropper. 

(v) In principle, one can also separate individual modes 
to match the requirements of the pulse user, e.g., to pre- 
pare a single mode maximally squeezed field, or to extract 
a maximally entangled two-mode field, etc. 

(vi) The pulse separation is generalization of spectral 
filtering which has also been used for observation of pho- 
ton number squeezing |l9[ . We have shown that spectral 
filtering can never beat the coherent amplitude modula- 
tion approach in achieving better photon number squeez- 
ing. In our approach, we can understand why spectral 
filtering can help with squeezing, but we can also see its 
limitations. 

(vii) Our approach of multimode description of solitons 
is different from the approach by Haus and Lai [Q (see 
App. [f] for more details). The formalism of JjJ was devel- 
oped to solve an idealized quantum soliton equation and 
to study the soliton dynamics. Its essential feature is the 
decomposition of the field into the "soliton" and "contin- 
uum" part, where the soliton field is completely described 
by four operators. Our approach ignores the dynamics 
(which is rather complicated in real life), focusing on the 
phenomenological description of the pulse. The four soli- 
ton operators of Q| can be expressed by means of our 
quadrature operators, but they alone appear to be in- 
sufficient for description of the observed phenomena. It 
is necessary to work with the complete set of "soliton" 
and "continuum" operators if one wishes to describe the 
observed pulses using the formalism of m. 
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APPENDIX A: SIMPLIFIED CALCULATION OF 
TRANSFORMED MEAN QUADRATURES AND 
VARIANCES FOR FEW EXCITED MODES 



APPENDIX B: QUADRATURE MOMENTS OF 
MULTIMODE GAUSSIAN STATES 



An JV-mode Gaussian state with mean quadratures \i 
and variance matrix V has the Wigner function W(ji), 



W(p): 



1 



(2w) N VdetV 



exp 



(Bl) 



Let us assume that only the first iV" modes of the b\ 
system are occupied, while the rest is in vacuum, \0) k for 
k > N. The quadrature means and variances are thus 



and 



Ji k = for k > N, 



V^k^k = ^ for k > N, 



(Al) 



Vnk,n'k' = for fi ± p! 

and (k > N or k' > N). (A2) 



The transformations (27) and ( |29| ) can then be written 
as 



and 



& = £<£*. 



(A3) 



where V~ x is the matrix inversion of V. Thus, a Gaus- 
sian state is fully determined by 2N 2 + 3N real parame- 
ters (2N mean quadratures plus N(2N + 1) independent 
elements of the symmetric matrix V). 

Generally, a quadrature moment (/}?) can be calcu- 
lated as the integral of the Wigner function 

</~4 l ) = J ... J tiWMdXx ...dP N . (B2) 

The symmetrical two- variable moments A ({Afc; Afc'}) and 
| {{ Afc > Afc' } ) can be calculated as 

1 



({Afc, A*/}) = 



H kl i' k ,W{n)dX x ...dP N (B3) 



and 



2 ({Afc, Afc'}) 



1 



.dP; 



A? 



6pk,n>k>, (B4) 



£ 



with 



4££^^> 4 ) 



M k>N 



Z 



A? 



E- E £ 

/i,fc n=X,Pk=l 



(A5) 



(A6) 



where d^n'W = if /Zfc and /ijy commute with each other, 
and S^k.fi'k' = 1 if Mfe and a4' are conjugate variables. Us- 



ing Eq. (Bl) one can analytically evaluate the integrals 
and express the quadrature moments by means of the 



parameters fi k an d V' 



We obtain 



fik,[i f k' ' 

(A*) = ~p k , 

(At) = vl + Vuk^k, 

(Al) = 7*fc + 3M fc ^fc,|ifc, 
(At) = /4 + 6 Afc^,M - 

({Afc, Afc'}) = AfeAfc' +^ufe,/u'fe'> 



3^ 



and 



(B5) 
(B6) 
(B7) 
(B8) 

(B9) 



To get (A5) from (A4), the orthogonality and complete- 
ness relation (20) was used. As a special case, let us 



consider the variance matrix elements VL A „ 4I when the 



mode functions /fc(w) are real, as in Sec. |IIB[ Eq. ( |A5| ) 
then becomes 



/ f » 2 - 12 1 \ -2-12 , - 2 t / i - '2 t / 

+V/ jl k,nkVf Jl 'k',fi'k' + ifikH'k'Vfik^'k' 



k,k'<N 



v'. , = ^1 



Note that although flB5|), (B6), and (p39|) are gen erally 
valid for all states, the equalities (B7), jBq), and (BK) 
only hold for Gaussian states. Corresponding expressions 
can be found also for the moments of the quadratures 
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APPENDIX C: PHOTON STATISTICS IN 
DIFFERENT FREQUENCY CHANNELS OF A 
SINGLE NONMONOCHROMATIC MODE 

Let us assume a nonmonochromatic mode defined by 
the discrete function f(uJk), k = 1 . . . N, J2k=i \f( UJ k )\ 2 — 
1. Let the state of this mode be Gaussian. We are inter- 
ested in the photon number correlations between different 
frequency channels. 

Theorem: All non-zero elements of a multi-mode pho- 
ton correlation matrix of an effectively single-mode Gaus- 
sian state have the same sign. They are negative (posi- 
tive) iff the single-mode state is sub(super)-Poissonian. 

Proof: The mode function f(u>k) defines the first row 
of the unitary transformation matrix U. Let the pa- 
rameters of the single-mode Gaussian state be X, P, 
Vxx, Vpp, and Vxp- The multi-mode parameters Xk, 
Pkj Vxk,xi, Vxk,ph etc., can b e cal culate d by means of the 
simplified summation as in (A3) and ( A5). The sign of 
the off diagonal element of the correlation matrix j36| ) 
is determined by the sign of the c ovar iance (|3^). After 
some algebra with applying Eqs. (A3), (A5), and (pi) 
we arrive at 



X 2 V xx - 



-\{Vxx-\ 



cov (n k ,ni) = \f(ui k )\ 2 \f(ui)\ 2 

2 / 

2" 



2XPV X p + P 2 Vpp- 



XP 



Uv FP -l 



.(CI) 



Thus, we can see that the sign of the non-zero elements 
does not depend on the arguments ojk and u>i , but is fully 
determined by the expression in the square brackets. If 
we denote the photon number in the single mode by h, we 
find that the quantity An 2 — (n) is equal to the expression 
in the square brackets. This quantity is negative for sub- 
Poissonian states and positive for super-Poissonian states 
by definition, QED. 



APPENDIX D: PROOF THAT THE MAGNITUDE 
OF THE SPECTRAL-FILTERING 
QUADRATURE VECTOR IS LESS THAN 1 



Showing that 

a 2 =y.{U 2 + pP)<i 



i<N 



(Dl) 



is equivalent to showing that 



^ CknQn (X k Xi + P k Pl) 

k,l,n<N 

< CkiiXkXi+PkPi), (D2) 



k,KN 



which follows directly from the definitions of X? and P? , 
Eqs. <M), @. To violate Eq. (EI), it would be neces- 



sary to have some vector Y such that 
Y T c 2 Y > Y T cY, 



(D3) 



where c is a matrix with the elements cu- Eq. ( p3| ) 
could only be valid if there exists some eigenvalue c\ of c 
such that c\ > 1. Let us assume that such an eigenvalue 
does exist. Let u k be the elements of the corresponding 
eigenvector, i.e., 



2J CklUi = C\U k , 
KN 



(D4) 



i.e. 



J^m; / c 2 {u)f k (uj)fi(uj)duj = c x u k . (D5) 



Let us define a function q(u>) as 



(D6) 



1<N 



Since Uk = J 9(w)/fc(w)(ia;, one has 

c 2 (u))q(u))f k {io)duj = c A / q{uj)f k {uj)duj, (D7) 



i.e. 



c 2 {Lo)q{u))fk{Lo)du> 



> 



q(u)f k {u))duj 



> (D8) 



which cannot happen for any function q(u>) since 
\c 2 {uo)\ < 1. Thus, A 2 < 1, QED. 



APPENDIX E: ELIMINATION OF COHERENT 
AMPLITUDES OF A MULTI-MODE GAUSSIAN 
STATE 

Let us assume that N modes of the 6-system are ex- 
cited and the rest is in vacuum. Then, one can redefine 
the N modes such that only one of them (say mode bi) 
has nonzero values of Xt, Pi, whereas X m = 0, P m = 
for m > 1 (the variance matrix elements corresponding 
to these modes are, however, generally non-zero). 

Proof: Define [3 m = 2~ 1 ' 2 {X m + iP m ) = (b m ), m 
— 1 . . . N. Let f m be the mode functions normalized 
such that (/ m ,/m') = 5mm 1 , where the bracket denotes 
the scalar product. Let us define a new mode func- 
tion gi as <7i = A^2^ =1 f3* n f m and a corresponding 
annihilation operator h\ = ^53m=i @mPm) where A — 

feLi IPml 2 )^ 2 - Then (hi) = Jj^LJ/U 2 - Let us 



define mode functions g 2 ■ ■ ■ gN as linear combinations of 
fi . . ./at by some orthogonalization procedure, i.e., g m 



= 2 m '=i Gmm'fm', m= 2...N such that (g m ,gm>) = 



5mm' : Tn = 1 . . . N. In particular, 



JV 

(5l ,9m) = Ay^fyGrrd 
1 = 1 



(El) 



13 



for m > 1. The annihilation operators corresponding to 

i'b m t, and their mean 

' T mm' 



G, 



the g-modes are h n 

values are (h m ) = ^\ n i—i G mm i (b m >) = 2~2m'=i G mm >f3ri 
= according to (JEl|) . Thus, in the new system of modes 
only the first one has a non-zero coherent amplitude, 
QED. 



APPENDIX F: RELATIONSHIP TO THE 
SOLITON PERTURBATION APPROACH BY 
HAUS AND LAI 

In |0] (see also |l^, [ll| ) it is suggested to describe the 
quantum fluctuations of the Nonlinear Schrodinger equa- 
tion (NSE) soliton by writing 



a(x, t) = a (x, t) + Ad{x, t), 



(Fl) 



where a (x) is the unperturbed solution of the NSE of 
the hyperbolic secant form, and 



Aa = Aa so i + Aa c 



(F2) 



is the quantum perturbation part with Aa so i describ- 
ing fluctuations of the soliton degrees of freedom, while 
Aficont corresponds to field fluctuations not contained in 
the soliton solution and thus belonging to the continuum. 
The soliton fluctuations are expressed by means of four 
operators An, AO, Ax, and Ap as 

Aa so i = [Aft(t)/„(s) + A0(t)f (x) 
( KA 2 \ 

+Ax(t)f x (x) + n Ap(t)f p (x)] x exp (t-g" 2 *) , ( F 3) 

where the functions /„, fg, f x , and f p are derivatives of 
the soliton function 



a (x, t) = A Q exp 



(KAl C , 

( -^t--plt+ Po x + e 



((h { ' .r ■'■^ , 



with respect to the parameters n (— 2£\A \ 2 ), p a , O , 
and x . In these equations the value n is the mean pho- 
ton number of the pulse, £ is the soliton width, p and 
x refer to the carrier frequency and soliton center po- 
sition, respectively. The value K is related to the Kerr 
nonlinearity and C to the medium dispersion. To find 
the values of the operators from measured data, one in- 
troduces adjoint functions / , / , / , and / with the 
property 



Re 



t k {x)fi{x)dx 



= S, 



k-l 



(F5) 



with k, I = n, 9, x,p (note that themselves do not form 
an orthogonal set; explicit form of functions /& and / 



can be found in p2|). Writing Aa as a sum of hermitian 
operators Aa = Aa^ + iAa^ one finds that 



An 
A9 
Ax 
Ap 



f n (x)Aa^(x)dx, 
i J Q(x)AaW(x)dx, 
J H(x)Aa^(x)dx, 
i— { f*(x)Aa^(x)dx. 



(F6) 

(F7) 

(F8) 

(F9) 
(F10) 



The relationship of this approach to our scheme can 
be easily examined if one assumes the first four mode 
functions /fc(w) in the form 



/iH = 
/ 2 M = 

/a(w) = 



sech — , 

'2uj LU 



3 u> . cu 

tanh — sech — 

2uj uj q lu 



(Fll) 
(F12) 



1 



1 



9 



. UJ u> u> , 
2 — tanh — sech sech 

UJ Q UJ a LU UJ, 



(F13) 



/ 4 (w) = i 



V3 



1 



tanh — sech — 

UJ Q UJ 



2 LO UJ 

sech — 

3 UJo UJ 



(F14) 



These functions were obtained by an orthogonalization 
procedure from the Fourier transformed functions f n ,e,x,p 
with uj — 2c/(7r£). With this choice of the mode func- 
tions one finds that the soliton perturbation operators 
can be expressed by means of the quadratures X\, X%, 
Pi, P 2 , P3, and Pi as 



An 

ao 

Ax 
Ap 



f 2n X l , 
1 



I2jn 
2c 



Pi + 



/ 2n n 



-Ph 



'§n uJo 



(F15) 
(F16) 
(F17) 



/ §n ~uJo 5 



2c 



y-^A. (F 18 ) 



Assuming vacuum fluctuations of X^ and Pk, {X 2 ) = 
1/2, (P 2 ) = 1/2, one finds that the uncertainty products 
of the soliton perturbation operators are 

(An 2 ) <A0 2 ) = I + ^Uo.675, (F19) 



3(i 



(Ax 2 )(n 2 Ap 2 ^ 



36 



0.274, 



(F20) 



14 



which are larger than the minimum uncertainty value 
1/4 following from the commutation relations [An, A9] = 
i and \Ax, n Ap] — i. This result corresponds exactly 
to that of [0, |l2). Our interpretation of the result is 
that the operator A8 is not purely a conjugate of An, 
but it contains an admixture o f an operator commuting 
with An. This admixture [in ( |F16| ) proportional to P 3 ] 
increases the noise above the minimum uncertainty limit. 



Similar interpretation holds for the pair Ax, n a Ap. 



From Eqs. ( F15 )-( F1S ) one sees that from the statis- 
tics of X± t 2 and Pi, 2, 3, 4 one can determine the statistics 
of the soliton perturbation operators An, AO, Ax and 
Ap. This does not hold vice versa: knowledge of the 
four soliton operators is not enough to determine the six 
quadratures. 
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